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The paper deals with the problem of surface effects at a fluid boundary 
q . produced by a step force field. A classical simple fluid with a locally placed 

• field simulating a solid is considered, 
i— — i" The specific surface 57-potential 7, the surface number density, and the 

Henry adsorption constant are determined. 
t-H \ A surface cluster expansion (the expansion of the surface omega poten- 
^ . tial in powers of the activity) is obtained. This expansion is similar to the 
■ cluster expansion for the pressure in which the integrals of the Ursell fac- 
t^J" ' tors are replaced by sums of the first-order moments of the Ursell factors 
ly-) ' over sectors or orthants. 

\Q [ The contact theorem is extended to the case of a finite step field. 

It is found that the surface number density (its invariant part) is de- 
1 termined by the first-order moment of the pair Ursell function taken over 
CO ■ the entire space. 

The high-and low-temperature limits are analyzed and are shown to be 
consistent with the previously obtained general results. 

Symmetry features of the solution with respect to permutation of the 
solid and fluid regions are established. 

X 

1. Introduction 



Statistical mechanics of fluids near an arbitrary permeable wall was constructed 
in pn. In the present paper, this theory is developed and applied to the adsorption 
problem with a finite step potential. In this case, an important role is played by the 
specially designed apparatus of sectorial and orthant moments. 

The importance of this problem is due to two things. First, the step potential is 
an important model system. Suffice it to say that this problem is an extension of the 
problem of a rigid impermeable wall and, in particular, it provides a generalization 
of the long and widely used contact theorem [2] (section [9]). 

In this model system, surface cluster expansion plays an important role, primarily 
as an accurate method for validating any results in the range of not too high number 
densities (sections [5l and appendix IB) . 
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Second, the step potential model is related to probabilistic problems. In particular, 
the behavior of the nonlinear surface coefficient at high temperatures is related to the 
surface part of the standard deviation of the number density (section [7]) . 

However, as noted in [3], even the general problem of a permeable wall had not 
been considered in terms of statistical mechanics until 1980. Later, the emphasis has 
shifted to approximate and computational methods. 

A number of results valid for the step potential is not considered in the present 
paper because they are formulated in general terms in HQ. These are, for example, a 
proof of the identity of the "tension" and "adsorption' approaches, a consideration 
of the "mechanical definition" of 7, and other issues. 



2. Primary relations 

2.1. Canonical ensemble 

The probability density of finding a given spatial configuration of a particular set 
of particles [H p. 181] is given by 

p [ k \ = 7T f exp(-PU 1 ... N )dr k+1 ...dr N . (1) 

V 

Here N is the number of particles in the system, f3 — 1/ksT, ks is the Boltz- 
mann constant, T is the temperature, E7j...jy is the particle interaction energy, and V 
the volume of the system. The integration is performed over the coordinates of the 
particles of the ensemble. Zn is the configuration integral 

Zn = J exp(-f3Ui... N )dri...dr N . (2) 
v 

Passing to the distribution functions for an arbitrary set of particles, we obtain 

(fc) _ Nl P w (3) 
Qc,i...k- (iv- /c)! 1 - fc " l j 

This quantity defines the probability density of finding a given configuration of k 
arbitrary particles for a canonical ensemble. 

2.2. Grand canonical ensemble (GCE) 

We average the equality over fluctuations in the number of particles, i. e., apply 
the operation X)tv=o ^° both sides of the equation. Here 

N <7 
y _ Z Z N 

Fn ~ N\E V (4) 

is the probability for the GCE to have a certain number of particles N within volume 
V. z is the activity: 

z = —fj-> ( 5 ) 
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where fj, is the chemical potential and A = h/y/2TtmkBT; h is the Planck constant, m 
is the particle mass, and Sy is the grand partition function of a system of volume V: 



— z"Z, 



- = i + E^ (•) 



N=l 



We obtain: 

oo 

= E *8?i...^. (7) 



JV=fc 

or 

vfe r °° „N 



(k) _ % 
0G,l...k ~ j; 



z N f 1 

exp (-/3Z7i...fc) + ^ — / exp(-j8t/i...jv+fc)drfe+i...drfc + jv • (8) 

AT 1 ' ** 



N=l y 



Qq , k dchne the probability density of finding a certain configuration of k arbitrary 
particles in the GCE. For an ideal gas, this quantity is equal to g k , where g — N/V 
is the number density. 



ik) 

The Ursell factors k are also called cluster functions. These functions appear 



2.3. Ursell factors 

The Ursell factoi 

in the well-known pressure expansion in powers of the activity (3J p. 135J3, [3 P- 230 
P(z,T) = zk B T + k B Tjr — y f U[ n) n dr 2 ...dr n n (9) 

n=2 H ' J 

They are also included in the corresponding expansion of the number density 

g(z) = z + zjr — I Ut + £ 1 dr 2 ...dr n+1 , (10) 
which is obvious from the formula 

-(f), 

taking into account the relationship ([5]). 

The Ursell factors have the locality property, which is important for our considera- 
tion: they decay rapidly as any group of particles, including unit groups, move away 
from the others. 

The Ursell factors can be defined by the equality [7] 



U^.k = Et- 1 )'"^ - X ) ! IT exp[-/3C/(n a )] 

w «=i (12) 

i 

1 < k a < k, ka — k, exp(—/3Ui) = 1, 



this paper, the expansion is not in integrals of the Ursell factors, but in their related Thiele semi- 
invariants. The relationship between them is defined by the transition from particular partitions 
to a partition topology (see below). 
2 The absence of the integration limits will always mean that the integration is over infinite space. 
The validity of infinite limits for this case is shown in \E\. 
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where {n} denotes some partition of a given set of k particles with coordinates r± , . . .r^ 
into disjoint groups n ai I is the number of groups in the particular partition, k a is the 
size of the group with the number a, and U(n a ) is the particle interaction energy in 
the group n a . The summation is over all possible partitions, and the meaning of the 
condition cxp(— (3Ui) = 1 is obvious: unit groups do not contribute to the product in 
this case. 
For example, 

1 

exp(-0C7i, 2 ) - 1 (13) 
exp(-/3C/ li2i3 ) - exp(-/3£/ lj2 ) - exp(-/3C/ li3 ) - exp(-/?*7 2 , 3 ) + 2. 



U 12 - 
M l,2,3 



2.4. Presence of an external field 

The configuration integral of a closed inhomogeneous system is given by 



N 

Z U N = J expi-P^m- pU 1 ... N )dr 1 ...dr N , (14) 

v i=1 

where Uj = u(rj) is the energy of interaction of the z-th particle with the field. 
For the GCE, we introduce the quantity 

JV=1 

which is obviously the grand partition function of the system in the presence of an 
external field. 

For arbitrary activities, the distribution functions have the form 



k K 

^... fe =ic/[n^]{ ex p(-^-* 



k k 

™%\\ exp(-/3U 1 ...k) 
L ;=i 

00 z N f r k+N i i 

+ E ]vT / i n exp ( ~~ /^i...A r +fc)dr fc+ i...(ir A r +fc j, (16) 



Af = l y i = k+l 

where 

f9i - exp(-0ui). (17) 
Varying the external potential, we obtain 

k 

S& = -P[&E Su * + I A<?Si)*«H-idTM-i]. (18) 

i=l ^ 

which for the case fc = 1 reduces to the Yvon equation 

Sqi = -P{qi5 Ui + J ^Su 2 dr 2 ), (19) 

A$ = Qil - Q1Q2 (20) 
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is a second-order Ursell function (in this case, for the system in the field ). 

If the variation of the potential is constant in spaceQ, it i s equal to the change in 
the chemical potential with the opposite sign, and from (|18|) we obtain 

/ - &il)dr k+l = z k+l i(fy\ (2D 

where the differentiation is performed at constant temperature. Thus, this relation 
has the same form as in the case of a homogeneous medium 0. 

Using the prop erties of the fractional generating function J3J) (see [I]), we obtain 
an analog of i|10|i 

oo n n+1 

e (n,z) = zOx + z6xJ2 Z — / [ II h]ut + nlidr2...dr n+1 , (22) 

where g{r\, z) is the number density in the presence of an external field. 
We will also need the following expansion: 

co n „ n+2 

jg>(*) = z 2 M 2 {z4 2) + £ Z - y \ [ I] e]u^ n %dr,...dr n+2 }. (23) 

n=l ' J 1=3 

An analog of (|23|) for a homogeneous medium is obtained in [5]. For an inhomo- 
geneous medium, it is sufficient to note that this relation is provided by the same 
recurrence relations for the partial localization factors pQ as in the homogeneous case. 



3. General case 



This section briefly summarizes the results obtained in [T] for the case of a permeable 
wall with the near-surface potential of arbitrary form. 
It is shown that 

Q u = -P( Z )V + [P(z) - P(z')} J ^(r)dr + v(z, z')A, (24) 

where O 17 is the omega potential of the system with the field placed in it, A is the 
bounding area of the field effect region, and 



v(z,z')= / e(x)P(z) + <p(x)P(z') -P*(x,z,z') dx 



(25) 



is a nonlinear surface coefficient (x is the coordinate directed along the field gradient). 
Here we have introduced the local pressure 



P*(r)= J Q(rW, 

— oo 

which, in this case, is equivalent to 



g(r) 



( dP*{r) 
\~diT~ 



(26) 



(27) 



3 It is sufficient that this region includes all the coordinates of the initial distribution function and 
that to their distance from the boundary of this region to the coordinates is at least a few atomic 
layers. The potential in this case may be arbitrarily inhomogeneous. 
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In addition, 

z' = zexp(-0u o ) (28) 

is the internal activity in the field region, u$ is the potential of the external field in 
the depth of the field effect region, and the quantities 

at \ _ exp [-/3u(s)] ~ ex P (~/M) , . 

K) ~ l-exp(-N {29) 

^ ~, v l-CXph^^)] 

1 - exp {-puo) 

generalize the Boltzmann factors (fTTjl and the Mayer functions 

tpi — 1 - exp (— jSttt) (31) 
respectively. It is easy to see that 

0; + ft = 1 (32) 

and that the functions 9i and ft decay in and out of the field, respectively, and tend 
to 1 otherwise. 

In addition, they have the property 

z9i + z'tfi — z6i. (33) 
The nonlinear surface coefficient is expanded in a scries in the activity as 

v{z, z') = k B Tj2-i 9^: n M n Lndr2...dr k , (34) 

n=2 ' J 

where 

+ r° 

dt'n = / {0(x)z n + if(x)z' n - [z6{x) + z'y(x)] J[ [z6(x + x t ) + z'£(x + a*)] }ds 

-co 1=2 

( 35 ) 

is a function of th e va riables X2, ■ ■ ■ ,x„ which is symmetric under permutations. 
The expression can be transformed to 

Q = -P(z)(V - V) - P(z')V + j(x')A, (36) 

where 

j(x') = - [P(z) - P(z')] (x' - x ) + v{z, z'), (37) 

x' is some arbitrary point within the transition zone or near it which defines the 
volume V , and 

oo 

xq = — J 8(x)dx + J (p(x)dx (38) 

-oo 

defines the zero a dsor ption plane (see below). 

The expression (l37)l gives the general form of the specific surface ^-potential. The 
last expression is divided into a term linear in pressure and the nonlinear term v. 

We see that there is a certain degree of arbitrariness in the division into volume 
and surface terms. However, for the step potential, its role is not as important, due 
to the presence of a distinct boundary. 
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Another form equivalent to (|37|) is 

x' -\-oo 

lW =/[P(/)-P-W]^ + /W-P-W] t (39) 

— OO x' 

Differentiating ()37|) with respect to the chemical potential and taking into account 
that the relation /rK<\ 

*-{£L (4o) 

retains its form, we obtain another expression for the surface number density for an 
arbitrary position of the separating surface: 

+00 

g s = [q{z) - g(z')] (x' -x )+ J [g{r) ~ 9{x)g{z) - y{x)g{z')]dx. (41) 

—00 

Here the parts linear and nonlinear in number density are clearly separated. When 
differentiating v(z, z') (and other quantities that depend on z') with respect t o z, the 
quantity z' must obviously be considered a function of 2, in accordance with (|28j) . 
Obviously, 

Qs = Qa,l + Qs,n, (42) 
where the linear part is defined by the equation 

Qs,t = [q{z) - g(z')] (x' - x ), (43) 

and the nonlinear part by 

+00 

Q a ,n= J [g(r)-6(x)g(z)-^(x)g(z')]dx. (44) 

— OO 

In addition, it is obvious that 

, J d[P(z)-P(z')] \ 
Qs.l = (x - x )^— ^ L j (45) 

and fdu(z,z')\ 

^ = -(-^1 • (46) 



T 



For the Henry adsorption constant 



K H = lim — (47) 
e^o g 

we have 

K H (x r ) = [x' - x (T)] [1 - exp(-^o)] (48) 

and the problem of determining Kh reduces to the calculation of xo . 

Note that Kh(xq) — 0, and the surface through xo can therefore be called the 
surface of zero adsorption. 

Using the expansion (|22[) . for the Henry absorption constant, we obviously have 

k H = lim ^ = exp(-/3u ). (49) 
e^o g(z) 
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4. Basic equations 



We now consider the step potential. 

We examine an open statistical system in a field such that the field effect area is 
smaller than the size of the system itself and the field is away from its boundaries; 
i.e., this is still the case of "a field in a system". 

The potential of the external fields is assumed to have a stepped form 

u(x) = { n° \ X< t Xstep \ (50) 

W \ (X> Xstep), y ' 

where x s t ep is the coordinate at which the potential energy of the particles has a 
step discontinuity, and the interparticle potential is arbitrary, provided that it decays 
rapidly enough for the convergence of the zero- and first-order moments of the Ursell 
factors. 

We introduce the characteristic functions 

1 (x < 0) 



= [ o (x > 0) ( 51 ) 

*<■)={; [iti <*> 

Then, for the step potential, we have 

6{x) = X(X ~ X s tep) (53) 

!p(x) = ip(x - x s tep) (54) 

and from (|25l) we obtain 

Xstep + 00 

v(z,z') = J [P(z') - P*(x,z,z')]dx+ J [P(z)-P*(x,z,z')]dx. (55) 



From (|38l . it is easy to find for this case that 

2^0 = Xstep, (56) 

and (j37f becomes view 

7 (z', z, z') = - [P(z) - P{z')] {x 1 - xstep) + »{z, z'). (57) 
The form of |39| . of course, remains unchanged. 

As in the general case PQ, for the step potential, we can make the substitutions 
P* —> P t and P* —> P s t in the basic equations |55)| and (f55)l , where Pt is the tangential 
component of the pressure tensor and P s t is the pressure on the transverse wall. 

For the surface number density, from ()41|l we obtain 

Xstep OO 

g s = [g(z) - q(z')] {x - Xstep) + J [g(x) - g(z')]dx + J [g(x) - g(z)]dx. (58) 

-OO Xstep 

In the case of the step potential considered in this paper, it is reasonable to choose 
the separating boundary x' to be equal to x s t ep or x s t ep ±D/2, where D is a parameter 
close to the particle diameter. The last variant is convenient because it includes the 
"dead" (conventionally speaking) volume, allowing the volume of one of the phases 
to be conserved. 
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For the Henry adsorption constant, from we have 

K H (x') = {x - x s tep) [l - exp(-/3u )] . (59) 

The meaning of the relation l|59|) is obvious: it is a coefficient proportional to the 
difference between the number densities in the limit of low activities [1 — exp(— /3uq)] 
multiplied by the volume of a parallelepiped of height [x' — x s t ep \ and unit area at the 
base. 

In this case, using (pf9"| , we arrive at the relation 

K H (x') = (x' -x st&p )(l-k H ), (60) 

between the Henry adsorption and absorption constants. Generally speaking, the 
relationship between these two constants is parametric - through temperature pQ , but 
in this case, it can be considered conventional. 

As in the general case of positive potentials [T], from (foTI)) it is evident that the 
adsorption and absorption in the range < ku < 1 are antagonistic. For a given x' , 
the adsorption is maximal at minimal absorption and tends to zero as ku — > 1 (100% 
solubility). 

For kn > 1 (negative potentials), the problem is redefined by symmetry (section HI 
and again reduces to the interval considered. 



5. Surface cluster expansion 

We calculate the n onlin ear surface coefficient using the first-order orthant moments 
technique (appendix IA.3|) . Our goal is to simplify the series ()34[l . or reduce it to a 
know n qu antity (this will be done in the next section). 

In l|34|l . passing to the function /( n_1 ) = / z n , we have 

y(z,z') =k B Tj2^ f A n -Ml. n dr2...dr n , (61) 

n=2 ' J 

where 



i=2 

Here 7 i 

\ = -= e - pU0 . (63) 
z 

Eliminating tp in view of ()32|) and using (|53|) . we have 

+°° n 

ft7n ] = I {A n + (l-A")x( a; )-[A+(l-A) x ( a; )]n[A+(l-A)x(a ; + a;j )]}^- (64) 

We expand the product and add and subtract x( x ) with the corresponding factor 
to each term, taking into account that 

+ oo k 

J \\\xi x + Xi) - x{x) dx = min(xi, . . . ,x k ). (65) 
»=i 
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The total coefficient of x( x ) vanishes and we obtain 

n 

=-J2 Xn ~ i ( 1 - X ) i E min^,^,...,^). (66) 

i— 1 samp 

The inner sum is taken over all subgroups of size i of the total group of n elements 
(zero plays the role of the first element). 

The Ursell factors are invariant under coordinate inversiorQ 

< = ~n, (67) 
where i = 1, . . . , n, which implies that for an arbitrary configuration, 

Z/W (r n ,r i2 ,..., r tn )=W« (-r n , -r h , . . . , -r in ). (68) 
Thus, we can equally well assume that 

n 

= E A " _l ( 1 - ^ E max(x Ql ,x Q2 , . ..,x ai ). (69) 

Denoting 



samp 



from (frJTjl . we have 



From (153, we obtain 



/n = E A "" i ( 1 - ^ \ I) ( J - !) / x 2 U [ ^ n dr 2 ...dr r , 



^ = J ft7nMi..ndT2...dr n , (70) 

oo n 

M*,A) = £V„. (71) 



j_2 ^ ' X2>0,Xs,...,X 

'n—\ 
i 



i I x 2 U { Q n l n dr 2 ...dr r . 



X2>X3,...,X i + 1 

+ (l-\) n (n-l) [ x 2 U$ n dr 2 ...dr n , (72) 

Jx2>0 7 X3,...,X n 

where the condition 

x ai ^> x a2 , . . . , X a j (^*^) 

implies the simultaneous satisfaction of the inequalities 

Xct\ ^ X&2 j • • • ) x ai ^> x a ^ . (*74) 

In d?2"jl , we used the fact that all groups of equal size give the same result and that 
unit groups vanish by antisymmetry. The presence or absence of zero in a group leads 
to terms of different forms, and the group of the maximum s ize always contains zero. 

Consider the first group of terms on the right-hand side of (f72"|) . Making the change 
of variables ^ r 2 (i — 2) 

r\-n (i = 3,...,n) (75) 



4 We use the quantum mechanical definition of the inversion [9] p. 96] 
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and taking into account that the Jacobian of this transformation is equal in modulus 
to 1 because of the triangular form of the functional matrix, we obtain (omitting the 
primes) 

/ x 2 U^l n dr 2 ...dr n = I x 2 U$ n dr 2 —dr n . (76) 

Jx2>0,x 3 ,...,Xi J x 2 >0,x 3 >0,...,Xi>0 

Here we us ed the invariance of the Ursell factors under the trans forma tion (fT5j) (see 
the property l|A.22|l for i =1, which is used in the transformation (|A.21|l ). 

The integral on the right-hand side of (|76| is not an orthant moment because the 
variables Xi+\ , . . . , x n may have arbitrary signs. But it is obviously the sum of orthant 
moments. That is, 

/ x 2^o n l..n dr 2-dr n = ( U ■ )[i+3,n]i- (77) 

J X2>0,X 3 ,...,Xi j_Q \ J / 

The second group of integrals on the right-hand side of (f72")l is divided into two 
parts x 2 > and x 2 < 0. With the first part, the problem reduces to the previous 
one (only the index i is increased by 1). 

The second part of integrals of this group (0 > x 2 > 23, . . . , £i+i) is transformed 
by the change of variables 



-r 2 (i = 2) 
Ti-r 2 (i = 3,...,n). 



(78) 



Performing similar calculations and using (|A.27|) for I = 2, we arrive at the expres- 



sion 

n—i—X 



L 



x 2 U$... n dr 2 ...dr n = - £ [ n \ 1 )[j + 2,n] 1 . (79) 

, 0>x 2 >x 3 ,...,Xi +1 j_q V J J 

Summation of (f72")l yields 

J n = (n-l)(l-A)^( n j(A"- fe - 2 -A fe+1 )[fc + 2,n] 1 . (80) 

k=0 ^ ' 

Multiplying |5D)l by z n , we see that this expression is symmetric with respect to 
permutation of z and z' 

z n I n = (n - l){z - z') y k J {z k+1 z ,n - k - 2 - z' k+1 z n - k - 2 )[k + 2, n] x . (81) 

Collecting terms with equal powers of A, we bring (|80[) to the form 
7„ = (n~l)(l-A){Kn] 1 (l-A"- 1 ) 

"-2 / n\ / o 

n — 2 \ , , , In — 2 



k=l 



+E( \ ) [B -*' B]l -U-i [ fc+1 '^) A 1- (82) 



The expression (fTTjl with /„ in the form of ((52"]) gives a variant of the surface cluster 
expansion. Other variants will be considered in the next section and appendix IBl 
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6. Surface quantities and the pair Ursell function 



f 2) 

Consider the first-order moment of J 7 ^ 2 (l20t 

m[ 2 \z,X) = /" x 2 T$dr 2 , (83) 

where the integration is performed over infinite space. Here Jy) 2 refers to the system 
in the step filed, and the x 2 axis is perpendicular to the field discontinuity surface on 
which th e co ordinate origin is located. 
Using (f2"3")l . we obtain 

oo „ 

Mi 2 \z,X) = e Y,T^\Jn, (84) 

/n 
x 2 [Y[e^ n l^dr 2 ...dr n . (85) 

We first calculate J„ - an analog of the J„ ((55) over the half-space x 2 > 0, 

/n 
x 2 [Y[9^ n l^dr 2 ...dr n . (86) 

Expressing the Boltzmann factors in terms of the characteristic functions (|51l) and 

Oi = tyi + Xi, (87) 
we obviously obtain 

n— 2 



where 



Returning to J n , we get 

J„ = J, i (A)-A"- 1 J„(A- 1 ), (89) 

where the second term on the right-hand side (the integral over the half-space x 2 < 0) 
is also determined by substituting ([57)) . 

Substituting ((55)) into (f59")) and comparing the result with ((52")) , we have 

/„ = (l-A)(n-l)J„, (90) 

or, in view of (|71l) and 



^(|^) =(l-A)flo 1 M«. (91) 

This is the relation of interest to us. It relates two fundamental quantities: the 
surface number density at the step boundary (see below) and the first-order moment 
of the pair Urse ll fu nctio n fo r this inhomogeneous medium. 

Substituting (|23)l into (|91)l and integrating over the activity at a constant temper- 
ature, we return to the expression of v as a series in z. Since the second-order U rsell 
function is determined by the number density and the pair distribution function (|20|) 
(i.e., quantities subject to constant attention ), this form is of particular importance 
in comparison with other variants (section [SJ appendix |B} . 
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The factor Oq 1 provides symmetry of the solution with respect to permutation of 
the solid and fluid regions (with a simultaneous change of the signs of uq and x and a 
change in the activity) , which must obviously be true for finite fields (see also section 
0). 

Indeed, if the observation point is to the right of the step (at the position +0), then 
#0 = 1 and we obtain 

=(l-A)M< 2) (z,A). (92) 



If the observation point is to the left of the step (set to —0), then #o = A and we 
arrive at the relation 

K£) = ~ x)M ' 2) ={1 ~ k )M< ? ){xz ' \ i (93) 

which is in agre eme nt with (|116|l . 

The equation (I91|l may be viewed as a surface analog of the compressibility equation 
for a homogeneous medium 



n^J, =M ° ' (94) 

Mg i) (z)=Jjf}dr 2 , (95) 

is the zero-or der moment of the pai r Ursell function. 
In view of the equation (|9Tj) can be written as 

fo>n = (A-l)^ 1 M 1 (2) . (96) 

The volume analog of ([9"6")) is transformed to 

e { e k B Tx T -l) = M (2) , (97) 

where kt is the is oth ermal compressibility. 

The expression ((9~lj) for A — > if the observation point is to the right of the bound- 
ary or for A — > oo if the obs erva tion point is to the left (low temperatures) is fully 
consistent with the results of |10] for a rigi d im permeable wall. This becomes clear if 
we take into account the form of the series (l23)l for this case. Furthermore, this result 
also corresponds to the surface parameter of an open stat istic al ensemble because of 
the identity of the probabilistic and potential constraints [TT] . 

7. High temperatures 

Consider the case of high temperatures (A — > 1). Let 

A = 1 + £ = 1 - pu + . . . , (98) 

then, from (|89|) . we obtain 



Jn - Ml) + eJ' n {l) - [1 + (n - l)e] [J„(l) - eJ^(l)] + . . . 

= £ [2j;(l)-(n-l)J„(l)] +..., (99) 

where th e prim e denotes differentiation with respect to the argument. 
From (]A.43[) it follows immediately that 

(n-2)J„(l) = 4^(l), (100) 
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and, thus, up to higher-order small quantities, 

Jn = ~nJ n (l), (101) 

or, in view of 

7„ = ^-n(n-l)J n (l). (102) 

Substituting (|102|l into (fTTjl and taking into account |2l{|) and (|56")l . we obtain the 
main result: 

v = ^ I x*4%dr 2 . (103) 

aj 2 >0 

This expression relates the nonlinear surface coefficient in the limit of small step 
size and the secon d-ord er Ursell function of a homogeneous medium. 

We show that (|103[) is consistent with the results of [T| for the potential of an 
arbitrary form. 

We start from the expression given in [T] , which in the case of a step potential takes 
the form 

|£) =(/3uo) 2 [ J x 2 F$dr 2 + J x 2 4% 3 dr 2 dr 3 ~\ + . . . (104) 

^ X2>0 X 2 >X 3 

Non-triviality is in the second integral on the right-hand side of (|104|) . which in- 
cludes the Ursell function of the third rank. To bring it to the desired form, we make 
use of the symmetry of this function with respect to coordinate inversion, permu- 
tations of particles, and translation. We consider two changes of variables with the 
transformation Jacobian of unit modulus. The first substitution, r' 2 — r 2 ,r' 3 = r 2 —r 3 , 
leads to the transformation 

J x 2 H%3 dr 2 dr a = J x 2 J r o 3 2 ) , 2 _ 3 rfr 2 o?r 3 

X2>X 3 X 3 >0 

= J x 2 F ( -l fi - 3 dr 2 dr 3 = J x 2 ^^ 3 dr 2 dr 3 = J x 2 T { ^ 3 dr 2 dr 3 . (105) 

x 3 >0 x 3 >0 x 3 >0 

The second change, r' 2 — r 2 — r 3 , r' 3 = — r 3 , gives 
J x 2 T ( ^ 3 dr 2 dr 3 = J (x 2 - x 3 ) F$_ 3 _ 3 dr 2 dr 3 

X2>X 3 X 2 >0 

= J (x 2 - x 3 )J 7 3 :3 l dr 2 dr 3 = J x 2 F$, 3 dr 2 dr 3 - J x 2 F^ 3 dr 2 dr 3 , (106) 

x 2 >Q x 2 >0 x 3 >0 

where in the last i ntegr al, we made one more change r' 2 — r 3 , r' 3 = r 2 . 
From (| 1 05|) and (|106|) , we obtain the relation of interest to us: 

X2>X 3 x 2 >0 

It has been shown previously [5] that 

j^i^-^U^f-)- (io8) 
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Thus, substituting (|107[) into (|104[) and using (|108|) . we obtain 



(dv\ (f3u ) 2 d f (2) 

Z 2 >0 

which is seen to be consistent with (|103|) . 

In 11 , it is shown that the first-order moment of the pair Ursell function defines 
the surface part of the rms fluctuation in the number of particles (at arbitrary tem- 
peratures). Thus, at high temperatures, the nonlinear surface coefficient is uniquely 
related to this quantity. 



8. Symmetry of the problem 

In contrast to the adsorption problem with an impermeable wall, the absorption 
problem, in principle, is symmetric. By varying the parameters, we may consider a 
particular region of the fluid to be free or situated in the field. This issue in a general 
setting is considered in pQ. 



Making the substitution 



z = z 
z = z 

x* = -x (110) 
u*(x*) = u(x) — Uq 
Uq = -u 

from ((2"9"|) and ((50"]) we make sure that 

6{u(x),u ] = 6[u*(x*) - Uq,-Uq] = ip[u*(x*),u^} , m , 
(p[u(x),u ] = !p[u*(x*)-Uq,-Uq] =6[u*(x*),Uq]. 

Thus, one can see, for exam ple, from (fM|) and (f55]l . the coefficient v is invariant 
under the transformation (|110|i 

v(9, tp\z,z') = v{tp* , 0*\z'*, z*) — v{9* , ip*\z* , z'*), (112) 
where we have introduced the notation 

£*=£[«*(**),«$] fn „s 
<y5 = <p[U [X ),U \, 

and the dependence of v on 9, (p in (|112|l is functional. 

In the derivation of (I112[) . it is necessary to use the invariance of the Ursell factors 
under inversion (|67|) . 

In the step potential problem considered in this paper, the specific form of the 
potential in the transition region gives rise to the additional condition 

0(x) = !p(-x), (114) 



which leads to a stronger consequence than (|112li : 

v(z,z') = v{z*,z'*). (115) 
This condition is equivalent to 

v(z,X)=v(z\,j) (116) 
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(cf. d)) or 

I n (X) = X n I n (h (117) 

(cf. d80|) V (As noted above, the equation (|116|l is consistent with the symmetry of 
the first-order moment of the Ursel l fun ction (section [6]).) 

In other words, the substitution (11101) in this case reduces to the symmetric form 

* / 

z = z 

z'* = z (118) 

z* = Xz 

X* = — (H9) 

* A 
x = —x. 



or the asymmetric form 

From (J55J it follows that 
and f5"7) implies that 



if we take into account that in the transformation 

X step = ~ x step- (122) 



q(x',z,X) = q(-x',Xz,j), (120) 
l{x',z,X) =j(-x',Xz, y), (121) 



The Henry constant (|4"5t satisfies the equation 

X 



K H {x',X) = XK H {-x',\) (123) 



in agreement with the general result pQ. 

9. Generalized contact theorem 

We employ the orthant moments technique (appendix [X| to generalize the contact 
theorem. Using this formalism, we first prove the usual contact theorem. 

Consider a distribut ion function of the first rank in the field of a rigid impermeable 
wall. As follows from (|22| . it has the form 



e (n) = z + V ^-j- / W 1 ( "+ 1 i +1 rfr 2 ...dr„ +1 , (124) 

n=1 n - Ja:2>0,...,x„ + i>0 

provided that the bounding field is in the range of negative x, and xi > 0. 

Let r± lie in the disco ntinuity plane of the potential; then in terms of the zero-order 
orthant moments (|A.9[l . this equation has the form 



z 

g(+0) = z + V— r [ n + l,n+l]o. (125) 

n=l 

Using the property (|A.15|) . we arrive at the series 



z 

? (+0)=^ + 2^^ T — [n+l] . (126) 

n— 1 ^ ' 
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Comparing this expression with (0, we see that the right-hand side of (|126|) is 
proportional to the pressure, and we obtain 

g(+0) = /3P(z,T), (127) 

the well-known contact theorem [12, p. 166]. 

We now turn to the generalized contact theorem. In the expression (f2"2"jl in this 
case, it is necessary to retain the Bol tzmann factors. Usin g the expression (l87)l and 
passing to the orthant moments (|A.9|) in ((2"2l . similarly to (|125|) . we obtain 

o(±0) ^ z n+l n / n \ 



6(±0) ^ n\ ^ \m 

n—l rn—0 v 

where the sign corresponds to the position of the observation point and where we 
have used the invariance of the Urscll functions with re spect to particle permutations 
and the symmetry of the zero-order orthant moments (|A. 12|) . 

Here 0(+O) = 1 

0(-O) = A 



(129) 



and, therefore, 

Using (|A.11|) . we arrive at 

e(±o) _ 



6(±0) 



+ E^r[ n + 1 < n + 1 ]°E xrn - ( 131 ) 

n— 1 m— 

Perfor ming th e inner summation and passing to integrals over the half-space by 
means of (|A.15[) . we obtain 

g(±0) z n+1 (l-A" +1 ) r . . 

Now we can again use the expression for the pressure (0; then, 

j^ ) =; + j r ^[i3P(i)-z-PP(^) + ^] (133) 

^-(AoW" (134) 

This expression is the desired result - the generalized contact theorem. 

In the limit A — > (low tem perat ures) or A — > oo (large negative potentials) , as o ne 
can easily see, the expression (I134|) transforms to the usual contact the orem (|127|) . 

For A — > l(the case of a ho mogeneous medium or high temperatures) ()134|) becomes 
the differential equation (fTTI) . 

It is easy to verify that the solution has the required symmetry property (|120|) 

q(±0,z,X) = q(tO,z\,j), (135) 

if we consider that 

0(-O,y) i 

l&ti) = A' (136) 
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Equation (|134|) under these assumptions can be written as 

g(+0) - e (-0) = P[P(+oc) - P(-oo)] , (137) 

where the pressure and number density are considered in the half-space of the corre- 
sponding sign. 

10. Summary 

1. For the step pote ntia l model, expressions were obtained for the specific sur- 
face H-potential 7<f5?| ■ the surface nonlinear coefficient ^ (|55|) . and the surface 
number density (|58|) . 

2. The zero adsorp tion surface in this case coincides with the discontinuity plane 
of the potential (f5tl|) . 

3. The Henry adsorptio n co nstant depends on the position of the separating bound- 
ary and is given by . 

4. The adsorption and ab sorp tion Henry constants are anticorrelated and their 
relationship is given by (|60[1 . 

5. Expressions were obtained for the nonlinear surface coefficient v as a se ri es in 
powers of the activity - an analog of the group series for pressure - (JZD, (f8"2l . 
section [B] and appendix [Bj The zero-order moments of the Ursell factors are 
replaced by sums of the first-order moments over orthants and sectors (appendix 

Ej. 

6. The invariant part of the surface number density at a stepped boundary is 
determined by the first-order moment of the pair Ursell function for an inho- 
mogeneous medium (section [B|. 

7. At high temperatures (small step size), the nonlinear surface coefficient v is de- 
termined by the first-order moment of the pair Ursell function of a homogeneous 
medium over a half-space (section [7]) ■ Thus, it is related to the surface part of 
the standard deviation of the number density. 

8. The step potential problem has the additional (compared to the general problem 
statement) property of symmetry with respect to permutation of the field and 
free fluid regions (section [5]). 

9. The well-known contact theorem is extended to the case of a step potential of 
arbitrary magnitude (section O. 

Appendix A Sectorial and orthant moments 

In this section, we consider the zero- and first-order moments of functions that have 
certain symmetry properties. Regardless of the specific form of these functions, it is 
possible to derive general relations which lead to a number of important results. As 
such functions we will use the Ursell factors. 

These moments will be taken over some regions: sectors and orthants (see below). 
The values of one integration variable (x) are restricted by certain conditions, and the 
other variables (y, z) are independent and take values corresponding to the infinite 
space. 

Sectorial moments play an auxiliary role and are required mainly to derive relations 
for the orthant moments. Ultimately, we are interested in integrals over whole space 
and half-space. 
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A.l Sectors and orthants 



The terms "sector" and "orthant" are used to refer to a consideration of a geometric 
picture in the space with coordinates X2, x%, . . . , x n . 
By a sector we mean a region 

x ai > x a2 > • • * ^> ^> ^> x a ^ 1 ^> • • • ^> x an: (A.l) 

where oti correspond to different particles and take values from 2 to n. 

The position of the zero in the chain (|A.1|) - j - can take values 1, . . . ,n for the 
zero-order moments and 2, . . . ,n for the first-order moments. It plays an important 
role and, in particular, determines the shape and size of the sector. The simplest 
example is given in Fig. [T] 



X 2 




Figure 1: Some sectors and the corresponding sectorial moments (of the first order) 
for n — 4. The moment-generating variable is xi- 

The orthant is the region 

x ai > 0, x a2 > 0, . . . , x aj _ 1 > 0, x aj+1 < 0, . . . , x an < 0. (A. 2) 

The total number of orthants is 2" _1 for the zero-order moments and 2" -2 for the 
first-order moments. 

Obviously, an orthant is the union of a certain number of sectors. The number of 
sectors in a given orthant is defined by permutations of particles within the positive 
and negative groups and it is equal to 

{j-l)\{n-j)\. (A.3) 

We now show that, due to the translational, permutation, and inversion symmetry 
of the Ursell factors, the sectorial and orthant moments satisfy some relations Q 



5 Of course, the relations obtained below are valid not only for the Ursell factors but also for arbitrary 
functions with the s ame symmetry properties. In particu lar, they are valid for the products of 
the Mayer functions Kit and for Ursell functions of type J20t . 
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A. 2 Zero-order moments 

Sectorial moments. In the zero-order case, the sectorial moments are defined by 
the equality 

(i)o= f U^ n dr 2 ...dr n , (A.4) 

J X2>X 3 >->Xi + 1 >0>X i + 2>->X n 

where i is the number of positive coordinates in the chain (|A. 1|) . (i = 0, . . . , n — 1). 
We note that the x coordinates of all particles should have a pronounced relation to 
zero. 

The change of variables 



(A.5) 



-r 2 (k = 2) 
rk~r 2 (k = 3,...,n) 

transforms the right-hand side of (|A.4[l as follows. 
The Ursell factor: 

W$,3,...,„ = (0, r 2 , r 3 , . . . , r n ) = «W (0, -r£, - r> 2 , ...,r' n - r' 2 ) 

= WW K, 0, r£, . . . , r'J = WW (0, r' 2 , r' 3 , . . . , <), (A.6) 

where in the penultimate equality we used the invariance of l/( n ) under translation, 
and in the last equality, its invariance under particle permutations. 
Sector (region of integration) : 

{x 2 > x 3 > ■ ■ ■ > x t+1 > > x l+2 > • ■ • > x n } 

= {~ x 1 > x 3 ~~ X 2 > ' ' ' x i+l — X 2 > > X i+2 ~ x 2 > ' ' ' > x n ~ x 2 } 

= {0>x' 3 >---> x' i+x > x' 2 > x' t+2 >■■■> x'J. (A.7) 

Given that the Jacobian of the transformation (1A.5[) is equal in modulus to 1, wc 
see that the zero-order sectorial moments are independent of the index i. Thus, we 
conclude that 

(*)o = (0)o, i = l,...,n-l. (A.8) 



m,n]o= / U { n l n dr 2 ...dr n , (A.9) 



Orthant moments. In the zero-order case, the orthant moments are defined by 

An) 

/, 0,2...n 

/x 2 >0,...,x m >0,a: m+ i<0,...,x„<0 

where the index m is equal to the number of positive coordinates plus one and n is 
obviously equal to the total number of coordinates (number of nonzero coordinates 
plus one), so that m = 1, . . . ,n. 

Obviously, an arbitrary orthant moment is equal to the sum of the corresponding 
sectorial moments: 

[m,n] = (m- l)!(n- m)!(0) , (A.10) 

which reflects the partition of the orthant into sectors (|A.3|l . This equality allows any 
moment to be expressed in terms of the higher-order moment 

n — 1 



[m,n] = [n,n] - ( A -H) 

\m — 1/ 

Note that (|A.ll|l implies the symmetry property 

[m, n]o = [n — m + 1, n]o- (A. 12) 
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For the integral over the half-space 

[n] = f U^l M dr 2 ...dr n , (A.13) 



X 2 >0 

we obviously have 

'n-2 

No 



m=2 



e(;: 2 )k-]o- (A.i4) 



From (|A.11[) . we obtain the relation 

r i x — ^ {rn — 1) r .. ?x r , / * \ 

No = 2^ ( x \ [ n Mo = -^[nMoi ( A - 15 ) 

m=2 ^ ' 

which, together with (jA.llj) . will be used in the derivation of the contact theorems. 

A. 3 First-order moments 

Sectorial moments. In the first-order case, they are defined by 

(i,n-j) 1 = I x l+l U { ^l n dr 2 ...dr n , (A. 16) 

J x 2 >->Xj>0>Xj + 1 >->x n 

where i is the positio n of the moment-generating variable (i = 1, ... , n — 1) in the 
chain of inequalities (|A.1[) and j is the position of the zero. Moreover, we assume 
that always i < j, i. e., the zero is always more distant in the chain than the moment- 
generating variable and that the latter is always positive. Thus, j = i + 1, ...,n. 
(The index n will always denote the rank of the Ursell factor.) Obviously, the second 
index, n — j, is equal to the number of negative coordinates. 

It follows that for an arbitrary sectorial moment of the first order 

(m,fc)i (A.17) 
the sum of the indices always satisfies the relation 

m + k<n-l. (A.18) 

Integrals of the type of (IA.16|) . where Xi+i < 0, are reduced to sectorial moments 
by the inversion ()67|l . 

Examples of the sectorial moments are 



or 



(1, 0)! = / x 2 l4j.. n dra ■ ■ ■ dr n (A.19) 

Jx 2 >->x n >0 

(2, n - 3)i = / x 2 l4%. n dr 2 ■ ■ ■ dr n . (A.20) 

The change of variables 



n+i (k = i + i) 

r i+1 -r k (k = 2,...,i,i + 2,...,n) 
transforms the right-hand side of (| A. 16|) as follows. 



(A.21) 
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The Ursell factor: 

= (0, -r[ + r' i+1 , . . . , r' i+1 , ...,-< + r' i+1 ) = U^{-r i+1 ,-r' 2 , . . . , 0, . . . , -r' n ) 

= U^(r l+1 ,r' 2 , ... ,0, ...,<) = WW(0, r£, . . . , r< +1 , . . . , r' n ), (A.22) 

where we have used the invariance oiU^ n ' with respect to translation, inversion, and 
permutations of particles. 

The sector (region of integration) : 

{x 2 > ■■■ > Xi> X i+ \ > X i+2 > ■ ■ ■ > Xj > > Xj+i > • • • > x n } = {—x' 2 + x' i+1 > 
• • • > — x i + x i+\ > x i+l > — x i+2+ a; i+l > ' ' ' > — ^j+'^i+l > > — a^j + i+^i+i > • • • 
> -Xn + x'i+l] = {-4 > • • • > -Xi > > -^ +2 > • • • > -Xj > ~x' i+1 > ~x' j+1 > 
■ ■ > -x'n) = Wn> ■■■> > > x 'j > > x 'i+2 > > x 'i > ■ ■ ■ > 4}- 

(A.23) 

Given that the Jacobian of the transformation (|A.21|) is equal in modulus to 1, we 
obtain the relation 

(*,n-j)i = (n-j + l,t-l)i, (A.24) 
where i — 1, . . . ,n — 1 ; j = i + 1, . . . , n or 

(m,k)i = {k + l,m- l)i, (A.25) 

where m = 1 , . . . , n — 1; fc = 0, . . . , n — m — 1. 

Equation (|A.25|) is the first basic relation for the sectorial moments of the first 
order. 

The change of variables 

P * = {r fc r -r, [k = l l... ,1-1,1 + 1,. ..,n), (A ' 26) 

where I is a fixed number in the range of j + 1, . . . , n yields a different transformation 
of the right-hand side of (|A.16|) . 
The Ursell factor remains invariant 

"o n 2 i n =U (n) (0,r 2 ,..., ri ,...,r n ) 

' = U^(0,r' 2 -r' l ,...,-r' l ,...,r' n -r' l )=U^(r' l ,r' 2 ,...,0,...,r' n ) 

= u^{Qy 2 ,...y l ,...y n ), (A.27) 

and the integration sector transforms as 

{x 2 > ■ ■ ■ > x i+ i > ■ ■ ■ > Xj > > Xj + i > ■ ■ ■ > xi > ■ ■ ■ > x n } = {x' 2 — x'i> ... 

> x' i+l — Xi > ■ ■ ■ > Xj — Xi > > Xj + i — Xi > ■ ■ ■ > —x[ > • • • > x n — x'{\ 
= {x' 2 > ■ ■ ■ > x' l+1 >--->x'j>x'i> x' j+1 >■■■> x'i_ x > > x' l+1 >■■■> x' n }. 

(A.28) 

In addition, in this case we obtain the transformation of the integrand factor 

Xj+i = x ' i+1 - x[ (A.29) 

and the integral on the right-hand side of (|A. 16|) splits into two. 

Given that the Jacobian of the transformation (|A.26[) is again equal in modulus to 
1, we obtain the second relation 

(i,n- j)i = (i,n- Z)i - (j,n- (A. 30) 
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where i = 1, . . . , n — 2; j = i + 1, . . . , n — 1; / = j + 1, . . . , n. 

Note that the index / is free, for given i, j, we have a whole class of relations (|A.30[) 
corresponding to different For example, for I = n, we obtain 

frn-jJi^vOji-ftO)!, (A.31) 

i=l,...,n-2; .7 = i + 1 n - 1. 

By sequential decomposition of (|A.30|) . it is easy to show tha t all t hese representa- 
tions are equiv alent and we can use any of them, for example, (|A.31|) . 

Thus, (|A.31I) is the second basic equ ation for the first -order sectorial moments. 

Substituting i = 1 and j = n — i in (|A.31|) and using (|A.25|) , we obtain 

(n-z,0)i + (z + l,0)i = (l,0)i, (A.32) 

where i = 1,2, . . . ,n — 2, which allows eliminating the higher-order terms (i, 0)i from 
the equations. Summing yields 

n — m 

£ (i,0)i = (|-m)(l > 0) 1 (A.33) 

i=m+l 

and, in particular, the total sum is 

n-i 

X)(<.0)i = ^(l,0)i. (A.34) 

i=l 

Orthant moments. In the hrst-order case, the orthant moments are defined by 
the equality 

[m,n]i= x 2 U^ n l n dr 2 ...dr n , (A. 35) 

J X2>0,...,X m ,>0,X m+ i<0, — ,X n <0 

where the index m is equal to the number of positive coordinates plus one and n is 
obviously equal to the total number of coordinates (the number of nonzero coordinates 
plus one). We will consider only the case where the moment-generating variable is 
greater than zer o, s o that m = 2, . . . , n. In the case x 2 < 0, we can use the inversion 
transformation (|67|1 . 

For example, the orthant moment 

[n,n]i= / x 2 U^ n 2 \ n dr 2 ...dr n (A. 36) 

J z 2 >0 : ...,2:, l >0 

specifies the first-order moment of the Ursell factor over the orthant where all Xi are 
positive. 

Obviously, the orthant and sectorial moments of the first order are linked by the 
relation 

m — l 

[m,n]i = (m — 2)!(n — m) ! (k,n — rn) i , (A. 37) 

fe=i 

which reflects a simple geometric partition of ortha nts int o sectors. 

Expanding the moment on the right-hand side of (1A.37|| with the use of (]A.3 1 f) , and 
performing the summation, we obtain 

m— 1 

[m,n]i = (m- 2)!(n - m)\ ^(fc,0)i - (m - l)!(n - m)!(m, 0)i, (A.38) 

fc=i 
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where m = 2, 3, . . . , n — 1. 

For m — n, we directly use (|A.37[) 



n-l 



[n,n]i = (n-2)!$^(M)i- (A.39) 
fc=i 

With the use of (|A.34I) . we make a further simplification 

^,^1 = ^^ (1,0)!. (A.40) 

Equation (|A.38|) defines an orthant moment in terms of sectorial moments. The 
inverse relation has the form 

( m 0)! = [rc-m + Mh K"]i + Mi (A 41) 

' n(m — l)!(n — m — 1)! n(m — 2)!(n — to)! n(n — 2)! 

where fc = 1 , . . . , n — 1 . The validity of (IA.41I) is easily proved by using the direct sub- 
stitution of (|A.38[) and the relation (| A.40|) . Thus, the sectorial and orthant moments 
of the first o rder a rc rigid ly con nected . 

By using (|A.38|) . (|A.40|) . and (|A.32|) and by direct substitution, we can verify the 
validity of the equality 



ra=1 



which will be employed to obtain high-tcmpcrature expansions. It will be convenient 
to bring to the form 

m=0 ^ ' m=l ^ 

by replacing ml = n — m. 



Appendix B Surface cluster expansion - third variant 

The ort han t mom ents are included in the surface cluster expansion for the step 
potential (f7T1) . I|82[l. Using the sorting identity (appendix [C|> . we can transfer the 
consideration to the level of sectorial moments. 

This identity (|C.6[) is applicable to the right-hand side of Using it, we obtain 

n 
i=l 

where the chain of elements is constructed so that 

xp t < xp s < ■ - - < xp n . (B.2) 

Naturally, the first element (zero) can be in an arbitrary position. Its position is a 
key factor. 

Since, a s al ready mentioned, the Ur sell factors are invariant under the inversion 
operation JSTJ, we can assume that [53] holds, or 

n 
i=l 



24 



and 

X 01 > X fo > ■ ■ ■ > X Pn- ( B - 4 ) 

Substituting the expression (|B.3|) into (|61|) . we obtain the quantity v of interest 
to us. In this case, since the integration is performed over the coordinates of each 
particle in an in finite space, it is necessary to consider all possible permutations in 
the chains (|B.4|) . bearing in mind that the zero can be in an arbitrary position. 

The importance of the position of the zero is due to the fact that it actually separates 
one chain into two independent chains. 

Obvio usly , for a given position of the zero, we have (n — 1)! identical integrals. 

Thus, ifCT)) becomes 

oo n n n „ 

pv{z, A) = (1 - A) J2 — E A 1 " 1 £ / X2 Utl. n dr 2 ...dr n . 

n=2 71 i=l I 1 * •■" '• 

,— I (B - 5) 

The second sum on the right-hand side of (|B.5|) is over the position of the integration 
variable in the chain of inequalities, and the third over the position of the zero, 
which vary from the first position to the last. The diagonal elements (in the case of 
coincidence of the positions of the zero and varia ble x%) are obviously zero. 

The integrals on the righ t-han d side of (|B.5[) . up to sign, are s ectori al mo ment 
of the first order (appendix IA.3|l . For i < j, we substitute (| A. 16|) into (|B.5|l . and 
for i > j we apply the inversion ()67[1 . Then, because in the t rans formation (I67L 
j' = n — j + 1; i = n — i + 1, we get for i > j that the integral in (|B.5|) takes the form 

- (n - i + l,j - l)i = -(j,n- (B.6) 

where the min us is due to the change in the sign of Xi . Here we have used l|A.25|l . 
Thus, from (|B.5|l . we obtain 

oo n n n 

Pv(z, A) = (1 - A) J2 Z ~ E E C * (B-7) 

n=2 i = l j = l 

where 

f (i,n-j)i (i<j) 
0^ = 1 (i=j) (B.8) 

I -{j,n- i)x (i > j)- 

Expanding the moments in (|B.8[l with the use of (|A.31[) and taking into account 
that [k, n— n)i = (&,0)i, we have 

f (i,0k <J = n) 

Cij = { (i,0)i-0',0)i (i,j^n) (B.9) 
I -(i,0)i (i = n). 



Summing yields 

oo n— 1 



/3v(z, A) = E z ™ [(! " A ) E Al_1 (*'' °)i - (! - A") 1 ^] , (B.10) 

n=2 i=l 

where we have used (|A.34|) . 

Taking out the factor (1 — A) from /„ |70| . we obtain 



(1 - A)n!{ E A* [H + 1, Q)i - - A"" 1 Mi} (B.ll) 

i=0 
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or, placing (1 — A) under the summation sign, 
/„ = „|{M>I + ]T A 1 [(i + 1, 0)! - (i, 0) x ] - A"- 1 ^ - 1, 0)! + A" -Mi}. (B.12) 

i=l 

Symmetry of this expression in the form ()117l) . as is easily verified, is provided by 
the relations (1A.32I). 

The expression (|B.11|) is equivalent to (f52")l . as is easily seen by equating the coeffi- 
cients of like p owers of A. Their coincidence is determined by the satisfaction of the 
relation (|A.41|) . 

Appendix C Sorting identity 

Consider the set of n elements to be compared (ordered) 

{xi,x 2 , ■ ■ ■ ,x n }. (C.l) 
Obviously, without loss of generality, we can assume that 

Xi < x 2 < ■ ■ ■ < x n . (C.2) 

We consider the structure 

n 

/ = ^(l-A) fe - 1 A"- fe+1 J2^M^ ai ,x a2 ,...,x ak ), (C.3) 

k=l samp 

where the inner sum is over all possible samples of k elements of n, and perform its 
calculation. 

Taking into account (|C.2|l . we have the expression 

n n— / ,\ 

/^(l-A^A"^ 1 E \T-i) Xl (C ' 4) 
fc=i i=i ^ ' 

and changing the order of summation, we obtain 

n n — l+1 / t\ n 

/ = r k Z J) (l-A)*^"" 1 = (C.5) 

1 = 1 k=l V" / l=\ 

or, ultimately, 

n n 

^A'x; = £(l-A) fc - 1 A n - fc+1 ^2mm(x ai ,x aa ,...,x a J (C.6) 

1=1 k=l samp 

if the condition (|C.2|) is satisfied. 

This is the desired identity, which may be called the sorting identity, because given 
the right-hand side of the equality, its left-hand side is automatically arranged in 
ascending order of xi . 

It generalizes the well-known maximum-minimums identity [13 a , which has the form 

n 

max(xi, . . . = y^(-l) fc-1 )] mm(x ai ,x a2 ,...,x ak ). (C.7) 

k= 1 samp 



Indeed, in the limit A — > oo, the equality (IC.6I) turns, obviously, in (|C.7|l 
As A — > or A — > 1, (|C.6I) becomes trivial relations. 
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Using the substitution x\ — —X{, we obtain the second form of the sorting identity 

n n 

^A'x / = ^(l-A)' £ - 1 A n - fc+1 ^x{x ai ,x a2 ,...,x ak ); 



1=1 k=l samp 

provided that: x\ > x-i > ■ ■ ■ > x n , 



(C.8) 



which can be called the inverse-sorting identity corresponding to the second equivalent 
form of the maximum-minimums identity 

n 

xmn(xi,...,x n ) =y^(-l) fe ~ x max(x ai ,x a2 ,...,x ah ). (C.9) 

k=l samp 
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